What happened to the gas-hquid transition in the system of dipolar hard spheres? 
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We explore the equilibrium properties of a system composed of dipolar hard spheres. A new 
theory based on the ideas derived from the work of Debye and Hiickel, Bjerrum, and Onsager is 
proposed to explain the absence of the anticipated critical point in this system. 



PACS numbers: 61.20.Gy; 75.50.Mm; 64.70.-p 

A representation of polar liquid in terms of dipolar 
hard spheres {DHS), N rigid spheres of diameter a and 
dipole moment p inside a uniform medium of dielectric 
constant constant eg is, probably, one of the most ba- 
sic statistical mechanical models. Yet, our understand- 
ing of this seemingly simple system is far from complete. 
A naive argument based partially on intuition and par- 
tially on oversimplified approximations suggests that as 
the temperature is lowered, a fluid composed of DHS will 
phase separate into a coexisting liquid and gas phases. 
This conclusion seems to be quite intuitive, after all if 
the potential between two dipoles, 

= , (1) 

is spherically averaged, one finds the familiar poten- 
tial of van der Waals which, of course, leads to phase 
separation. This argument, however, does not withstand 
the test of computer simulations which, up to now, have 
failed to locate any vestige of phase transition (2|-U] . In- 
stead, the simulations find that, as the temperature is 
lowered, the dipolar spheres associate forming polymer- 
like chains Can the formation of chains explain the 
disappearance of the liquid-gas transition? 

To respond to this question is not easy. In search for 
the answer it is interesting to recall the mechanism of 
phase separation in a different, but very related system, 
the restricted primitive model (RPM) of electrolyte . 
In that case, ions are idealized as hard spheres half of 
which carry a positive charge, while the other half carry 
a negative charge. At low temperatures, formation of 
clusters composed of positive and negative ions is ener- 
getically favored. First, appear dimers made of -I — pairs, 

then trimers -I h, etc. This looks very similar to the 

formation of chains in DHS, and yet the RPM does phase 
separate, while the DHS do not. What is responsible for 
this fundamental difference? At first, one might try to 
appeal to purely electrostatic considerations. Thus, it is 
tempting to attribute the phase transition in the RPM to 
the fact that by the time a cluster grows to contain four 
ions, the linear configuration becomes energetically unfa- 
vorable, and ions tend to arrange themselves in a square. 



These compact configurations could, in principle, provide 
the nuclei for the start of condensation. It is tantalizing 
to think that this is the essential difference between the 
ions and the dipoles; ions energetically prefer compact 
clusters, while dipoles prefer linear chains. As appealing 
as this argument might sound it is, nevertheless, incor- 
rect. A careful analyses of energies clearly demonstrates 
that the compact configurations also become energeti- 
cally favored for DHS by the time the clusters grow to 
contain four or more dipoles §. Thus, a square cluster 
in which the dipoles are arranged circularly (head-to-tail) 
at forty five degrees to the lines connecting the centers 
of the nearest neighbors, has electrostatic energy smaller 
than a linear chain of dipoles. The fact that the simu- 
lations observe polymer-like chains of dipoles, instead of 
compact configurations favored by electrostatics, implies 
that the entropy plays an essential role in formation of 
clusters. The energetics alone is insufficient to explain 
the distinct thermodynamic behavior exhibited by the 
RPM and the DHS. 

The analogy between the RPM and the DHS sug- 
gests that the methods developed to study the RPM 
might also be applicable to the exploration of DHS. In 
this respect the Debye-Hiickel Bjerrum theory DHBj 
of electrolyte has proven particularly illuminating [^,^. 
This theory augments the idea of screening, introduced 
by Debye and Hiickel , to explicitly take into account 
the formation of clusters composed of oppositely charged 
ions The estimates of critical parameters based on 
DHBj theory are, so far, the closest to Monte Carlo sim- 
ulations While the idea of cluster formation is di- 
rectly applicable to DHS and has already been exploited 
by various authors dipoles, unlike ions, do not 
produce any screening. Instead, the thermodynamic ef- 
fect of dipolar motion translates into renormalization of 
the effective dielectric constant of the medium. The ques- 
tion that we would like to answer is whether this residual 
interaction is sufficient to produce phase separation. 

We shall proceed in the spirit of DHBj theory ||] . The 
reduced free energy density, / = (SF/V , of solution will 
be constructed as a sum of terms embodying the most 
relevant physical features of the system, starting with 
the entropic ideal gas contribution /*'^ = p\xi{pK^) — p. 
Here p = N/V is the density of dipoles, (3 — l/ksT, 
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and A is the thermal wavelength. To obtain the electro- 
static free energy let us fix one particle at the origin and 
calculate the electrostatic potential that this dipole feels 
due to the presence of other molecules. The electrostatic 
potential can be found from the solution of Laplace equa- 
tion, = 0, supplemented by the appropriate bound- 
ary conditions. We shall separate this potential into two 
parts, (jjin for r < a, and c/fout for r > a. Clearly r < a 
corresponds to the excluded volume region into which, 
due to the hardcore repulsion, no other particles can 
penetrate. The boundary conditions require continuity 
of the potential, 0i„(a) — 4>out(o), and the displacement 
field, eo(/)-„(a) = e^o„j(a), across the surface r = a. We 
have introduced the renormalized dielectric constant of 
the bulk e, the expression for which can be obtained from 
the Onsager's reaction field theory p7[, 



(e-eo)(eo + 2e) 



(2) 



The Laplace equation can now be integrated to yield the 
potential of the central dipole due to other particles. The 
electrostatic free energy of the whole system is obtained 
through the Debye charging process [p|pO|] in which all 
the particles in the system are charged simultaneously 
from zero to their final dipolar strength. 
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dd 
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(3) 



The integration can be done explicitly yielding the elec- 
trostatic free energy density Q, 



- In f —-4 ) + 3 In V'(m) 



with 



e(u)/eo EE = i (1 + u) + ^^9 + 2^ + 1*2 , 



(4) 



(5) 



and u = 47rp*/r*, where we have introduced the re- 
duced density p* — pa? and the reduced temperature 
T* = ksTeoa^ /p'^. Combining Eq. (||) with the entropic 
contribution mentioned earlier, the total free energy den- 
sity of the system becomes / — /"^ -t- f^"^. It is a sim- 
ple matter to see that as the temperature is lowered, 
this free energy violates the thermodynamic convexity 
requirement, what results in a phase separation into co- 
existing high and low density phases ||l^,|l9| . Specifically, 
we find the critical parameters to be, p* — 0.0390807 
and T* — 0.138904. In principle, we could have also in- 
cluded the excluded volume contribution to the total free 
energy, expressed through the free volume or Carnahan- 
Starling approximation, but this would not significantly 
affect the location of the critical point |2^. The funda- 
mental conclusion of this Debye-Hiickel-Onsager theory 



[DHO) is that the system of dipolar hard spheres sep- 
arates into a coexisting liquid and gas phases. Can this 
result be trusted? Clearly, based on our experience with 
the RPM [H this conclusion must be taken with a grain 
of salt. Just like pure DH, the DHO theory is linear. 
This means that, although the DHO is quite adequate for 
capturing physics of large length scales, it fails for short 
distances. In particular, the DHO theory does not take 
into account the low temperature propensity to form clus- 
ters. It is precisely the importance of these configurations 
which is lost in the process of linearizations [ p^ leading 
to the Onsager relation (|^). This conclusion is very sim- 
ilar to the one reached for RPM |^. A solution, in that 
case, has been proposed more than seventy years ago by 
N. Bjcrrum, who suggested that the non-linearities, in 
the form of clusters, can be reintroduced into the DH 
theory through the allowance of "chemical" association 
between particles ]ll[ |. A theory based on Bjerrum's con- 
cept of chemical equilibrium has proven quite successful 
at treating the phase separation in RPM This sug- 
gests that the same kind of methodology might also be be 
useful for studying DHS. We, thus, suppose that at low 
temperatures the system consists of some free unassoci- 
ated dipoles of density pi , as well as clusters containing 
2 < n < oo hard spheres. The density of a n-cluster is 
Pn- The particle conservation requires that 



, npr, 



(6) 



Following Bjerrum, we shall first treat clusters as non- 
interacting ideal species. The interactions, therefore, are 
restricted to unassociated dipoles and their contribution 
to the total free energy density is given by Eqs. (|) and 
(|), with u = 4npl/T*. In the case of RPM model, this 
approximation has proven to be sufficient to locate the 
critical point |^ . The free energy density of a n-cluster 
reduces to the ideal gas form. 



/;'' = p„ln(p„A^"/Cn)-A 



(7) 



where we have introduced the internal partition function 
of a n-cluster. 



(8) 



Here C/f''* is the pairwise interaction potential obtained 
from p), and f2's are the relative angular orientations of 
dipoles forming a cluster. In the limit of low tempera- 
tures, where the DHO predicts the location of the critical 
point, the integrals in ( p|) c an be evaluated for chain- like 
configurations to yield p^j9[|. 
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where ijj^^^n) and tp^'^'>{n) are the polygamma functions 
of the first and second order, respectively. The condition 
for chemical equilibrium between dipoles and clusters is 
expressed through the law of mass action, /i„ = rifii, 
where the reduced chemical potential of a species s is, 
/is = df/dps- Substituting the total free energy density, 
/ — J2^=i frf + f'^'^^ fiii'i tlis distribution of cluster 
densities to be 



/- n nn^ 
Pn = CnPl e '^1 



(10) 



where the excess chemical potential is defined in terms 
of the excess over the ideal gas contribution, in this case, 
f^T — df'^'^/dps- It is important to note that within the 
Bjerrum approximation the excess chemical potential de- 
pends only on the density of free dipoles, and the expres- 
sion ( [lO| ) reduces to an infinite set of decoupled algebraic 
equations. 

We now make the following fundamental observation; 
since the clusters are ideal, their presence can only shift 
the critical density, while leaving the critical tempera- 
ture unaffected [^. Thus, the critical point must still be 
located at T* — 0.138904 and must still have the den- 
sity of free dipoles pl^ = 0.0390807! The distribution 
of clusters at criticality is obtained by substituting these 
parameters into Eq. ([l0|). In order for the sum in (^) to 
converge, the Cauchy-Hadamard theorem requires that 
A = lim„^oo pl/^ < 1. Inserting the critical parameters 



into Eq. (|10[), we find that at criticality Ac « 100 and 
the theorem is strongly violated. The critical density p^^ 
lies far outside the radius of convergence of (|). This 
means that for any finite density p, the density of free 
dipoles remains insufficiently small to reach phase sepa- 
ration! Clearly, the argument presented above assumes 
that only free particles interact while the clusters are 
treated as non-interacting ideal species. This certainly is 
a very strong approximation which must be considered 
in more detail, nevertheless, we note that a similar ar- 
gument has proven to be sufficient to locate the critical 
point of the RPM ||^. In that case, it was found that in 
the vicinity of the critical point, the series (H) was very 
quickly convergent with most of the ions belonging to 
dipolar pairs 

To explore the role played by dipole-cluster and 
cluster-cluster interactions it is necessary to account for 
their contribution to the overall free energy. This is far 
from simple. Some progress, however, can be made if 
we make the following observation. The electrostatic po- 
tential produced by a rigid line of dipolar density p/a 
is exactly the same as the potential due to two ficti- 
tious monopoles of charge ±p/a located at the line's ex- 
tremities. This can be shown explicitly by integrating 
Eq. (|l|). The isomorphism between line of dipoles and 
two discrete monopoles suggests that for low tempera- 
tures, when the dipolar chains are quite rigid, the dipole- 
cluster and cluster-cluster contribution to the total free 
energy can be approximated by the energy that is re- 
quired to solvate Nc — 2V J2^=2 Pn monopoles in the sea 



of dipoles, and by the energy of their mutual interaction. 
The solvation energy of an ion can be obtained follow- 
ing the same method presented earlier for calculating the 
dipole-dipole contribution. We find. 
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Performing the integration, the reduced free energy den- 
sity due to dipole-cluster interactions is found to be 
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(12) 



Finally, the cluster-cluster contribution can now be es- 
timated as the energy of a plasma composed of Nc ions 
inside a medium of dielectric constant e. We find the 
familiar Debye-Hiickel expression ||,|l^ , 
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where now 
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It is easy to check that at low temperatures both dipole- 
cluster and cluster-cluster contributions are quite small, 
and are unlikely to modify the previous conclusion of the 
absence of criticality in DHS. The exact calculation is 
rather difficult to perform since the law of mass action, 
when the dipole-cluster and cluster-cluster interactions 
are included into the total free energy, reduces to an infi- 
nite set of coupled algebraic equations. The preliminary 
analysis of these, based on a variational approximation 
for the distribution of clusters, does not, however, find 
any indication of phase separation. The details of these 
calculations will be presented elsewhere. 

We conclude that the low temperature propensity to 
form weakly interacting clusters absorbs most of the 
dipoles, preventing the density of free unassociated par- 
ticles from reaching the minimum necessary for phase 
separation. 
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